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We review recent developments in lunctional renormalization group (RG) methods for interact- 
ing fermions. These approaches aim at obtaining an unbiased picture of competing Fermi liquid 
instabilities in the low-dimensional models like the two-dimensional Hubbard model. We discuss 
how these instabilities can be approached from various sides and how the fermionic RG flow can be 
continued into phases with broken symmetry. 

I. INTRODUCTION 

Competing Fermi liquid instabilities are a general feature of low-dimensional interacting lattice electron systems. 
These instabilities are driven by long-lived particle-particle or particle-hole excitations near the Fermi surface (FS). 
They can lead to long range correlations of various types and a destruction of Fermi liquid behavior. In two and 
higher dimensions perturbative renormalization group methods are probably the least biased approach to study such 
systems at weak coupling and low energy scales in a transparent way. 

In the last decade various functional renormalization group (fRG) methods Q, S H S for many-fermion sys- 
tems have been developed. These schemes are capable of capturing the interplay of various low-energy scattering 
processes within certain approximations. However, the practical experience shows that each scheme has advantages 
and disadvantages. This has led to the development of alternative flow schemes|^ Q- They provide the possibility 
to tackle the competing instabilities from different angles. Fortunately, it turns out that besides well-understood 
conceptual differences the flow of these schemes agrees quite favorably over a wider parameter range of models like 
the two-dimensional Hubbard model. With this understanding we hope that further improvements of the approach 
are on a firm basis. 

This short review is organized as follows. First we outline the general formalism of the renormalization group 
for the one-particle irreducible vertex functions]^ 0, 0]. We explain how cutoff flow, temperature flow (T-flow) and 
interaction flow (IF) are obtained. Then we review some results for the two-dimensional Hubbard model. Finally we 
describe a novel approach how the fermionic RG flow can be continued into phases with broken symmetry. 

II. FUNCTIONAL RG FOR IPI VERTEX FUNCTIONS 
A. General setup 

The class of models we study here is described by an action of the form 

+ JN ^(f'l ^ + q)'ipk+q,s'>Pk'-q,s'tpk' ,s''ipk,s ■ ^^^^ 

Here i^{k, s) and ip{k, s) are Grassmann fields representing fermions with wavevector fc, Matsubara frequency fco (we 
write k ~ {k(),k)) and spin projection s = ±1/2. Q{k) denotes the quadratic part of the action, here given by 



Q{k)^T -ika + e{k) 



For fermions on a 2D square lattice with nearest neighbor hopping t, next-nearest neighbor 



hopping t' and chemical potential fi, the dispersion is 

e(fc) — ~2t (cos kx -f cos ky) — At' cos k^ cos ky — ^ . (2) 

V{k, k' , q) defines the spin-rotationally invariant, frequency- and wavevector-conserving interaction between two 
fermions. For the repulsive onsite interaction of the Hubbard model V{k, k', q) = U > 0. 

The connected correlation functions of the theory defined by Q can be obtained by taking derivatives of the 
generating functional, 

e-wiu) ^ f s)V,pik, s) e-^('^''^)+S'..= [^(k,.)ak,s)+ak,s)i,ik,.)] 



FIG. 1: RG equations for the two-point and the four-point vertex. The slashed line denotes a single-scale propagator Se{p). 
The one-loop graph for the four-point vertex 7^'*-' includes particle-particle and particle-hole contributions. With our truncation 
7^™' — for m > 6 the feedback of the 7^^' on 7^*' is neglected. 



For our purposes it is simpler to work with vertex functions, which are generated by the Legendre transform of W{^, ^) 
with respect to (p — SW/S£^, 



r(0, $) = Wit - E ['^(^' ^)^(^' ^) + ^"(^' *)'^(^' ^)] 



k.s 



This functional can be expanded in monomials of its sources, 



m! 

m>0 K 



(4) 



(5) 



_K is a multi-index which contains m frequencies, wavevectors, spins and Nambu indices of the $-fields. Nambu index 
-|- stands for and — for cj). The 7^™) are the m-point vertex functions of the theory. They are antisymmetric with 
respect to interchange of two particle coordinates. 

The two-point vertex 7^^-' is the inverse of the full propagator. For spin-rotational and U(l) invariance the antisym- 
metric four-point vertex for incoming particles fci, si, ^2, S2 and outgoing particles fca, S3, k^, S4 (the quantum numbers 
of particle 4 are dictated by the conservation laws) can be expressed in terms of a coupling function V{ki, fc2, fcs) byQ 



'S1S2S3S4 



{ki,k2,kz,ki) 



V{ki,k2,k3)6s^s3^s2Si 



V{k2,ki,k3)Ss^sJs 



(6) 



Now let us assume that the quadratic part of the action depends continuously on a parameter £ while the interaction 
terms does not contain £. Then we can derive an exact equation 5, 8, 9] for the change of Tiicj), 0) when £ is varied by 
taking the ^-derivative of Eq. 0] Inserting the field expansion (jSJ yields an infinite hierarchy of differential equations. 
This set of equations defines the renormalization group flow of the IPX vertex function with respect to the parameter 
£. In order to make this system tractable, one has to truncate it at some point. Here all IPI vertex functions with 
TO > 4 are set to zero. This is a first approximation that is common to all RG schemes discussed below. The remaining 
equations for 7|^-' and ji'^hl are shown graphically in Fig. 

We emphasize that in the exact hierarchy the right hand side contains only one- loop diagrams. This is one of the 
blessings of the approach which in certain cases even allows for a non-perturbative treatment 13] . Higher- loop diagrams 
are included via the ^-dependence of the vertices. They can be made explicit by reinserting the integrated flow of the 
vertices into the right hand side'lO'|. The equation for the coupling function Vi{pi,p2,P3) with 734 = pi + P2 — Ps isQ 



-^Veipi,p2,P3) = Tppj + 7^^, 



PH.l 



(7) 



with the one- loop particle-particle contributions Tpp £ and the two different particle- hole channels Tp^^ ^ and i 
where 



Tpp^i{pi,P2;p3,P4:) = 

- J dpVt{pi,p2,p) L{p, -p + pi + P2) Vi{p, -p + pi +P2,P3) 
TphApi^P2]P3,Pa) = 

dp -2Ve{pi,p,p3)L{p,p + pi -P3) Ve{p + pi -p3,P2,p) 
+Ve{pi,p,p + pi -p3)L{p,p + pi -p3)Ve{p + pi -p3,P2,p) 
+ Ve{pi,p,P3) L{p,p + pi -p3)Vt(p2,P + Pl -P3,p) 



(8) 



(9) 
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dpVi{pi,p + P2 -p3,p)L{p,p + P2 - P3)Vi{p,P2,P3) (10) 

In these equations, the product of the two internal lines in the one-loop diagrams is 

L{p,p') = S,{p)W^^\p') + wPip)S,ip') (11) 
with the so-called single-scale propagator 

Si{p) = -WP{p) 



WPip). (12) 



(2) 

Here, {p) denotes the full Green's function at RG scale 



B. Momentum-shell RG, temperature flow and interaction flow 

So far the setup of the fRG equations has been very general. Now we specify the flow parameter We discuss three 
different choices (although many more are possible). The guiding principle for a good flow parameter is that it should 
enable us to approach a specific singularity in the perturbation expansion in a controlled way and h) to include 
all other (possibly also singular) tendencies during the flow. Then the RG will give a more realistic picture than 
perturbation expansions which single out one dominant channel. In bubble or ladder summations, the singularities 
normally arise due to a pile-up of logarithms. These are roughly given by g log[M^/max(T, A)] with a coupling constant 
f/, bandwidth W and lower energy cutoff A. Hence, if we want to build in such dangerous terms step by step, we can 
either vary A, the temperature T, or the coupling g. This is made more precise below. 

a. Cutoff-RG Here we introduce a cutoff-function CA[e(fc)] into the quadratic part of the action, 

QA(fc) = T ^ CMk)] ^fc [-^fco + e(fc)] ^fe ■ (13) 

CA[e(fc)] is very large for |e(fc)| < A and CA[e(fc)] = 1 for |e(fc)| > A such that modes below A are not integrated over 
in the functional integral. In practice one mainly needs the inverse CA^[e(fc], which can be chosen conveniently as a 
sharp step function for analytical manipulations or as a smoothed step function for numerical treatments. The full 
Green's function is suppressed for modes with |e(fc)| < A. A frequency cutoff (or combinations of frequency and band 
energy) falls into the same class of cut-off RG. For our fc-space resolved approach described below a pure energy-cutoff 
is however most convenient. The momentum- shell RG is the widely used standard^ for many-fermion systems. It 
can also be derived from other exact RG equations [3. 1^ H. ITllIT^. It approaches Fermi liquid instabilities for fixed 
model parameters and fixed temperature. One of the successes of the method in the weakly coupled 2D Hubbard 
model are the clear signatures for d-wave superconductivity over a wide parameter range fS. l4Lll6l[. 

A serious drawback of the cutoff RG is however the non- uniformity in the RG scale at which one-loop particle-hole 
(PH) processes with different wavevector transfers q are included at low T. Typically PH pairs with a particle k and 
a hole k + q and large \q] ^ kp are integrated out at all scales A in the flow, depending on how close k and fc -I- ^ 
are to the Fermi surface. In contrast with this, PH processes with small q can only occur in the vicinity of the Fermi 
surface. For T the support of the PH bubble for w = 0, g — * shrinks to a temperature-smeared 5-function on 
the FS with width ~ T. Thus, the q* — > PH-modes are integrated out only for cutoff A < T, even if the density 
of states is divergent, and these processes give a singular contribution for T ^ 0. However, in a coupled flow with 
a lot of other tendencies, e.g. the omni-presence of Kohn-Luttinger superconducting instabilities, the flow normally 
diverges before we get down to A ~ T. The flow has to be stopped and the (f — *■ PH pairs did not have the chance 
to contribute to the flow by construction. Therefore, the cutoff-flow is still biased. Other approaches are needed to 
study the influence of the q particle-hole excitations. 

The following two schemes work without a cutoff and allow for a uniform inclusion of the various one-loop processes. 

h. Temperature flow The temperature-flow scheme^ uses the temperature T appearing in front of the frequency 
sum and in the Matsubara frequencies as flow parameter. Before the T-derivative can be taken, the fermionic fields 
need to be rescaled in order to remove the term appearing in front of the interaction term in the original action. 
Then the strategy is as follows. The theory is defined at some high T ~ Tq in terms of its two-point and four-point 
vertex. Here, we simply choose the free propagator and the local Hubbard repulsion as the initial values at Tq. This 
should be reasonable for sufficiently high Tq ~ bandwidth, as all perturbative corrections decay with a negative power 
of T. In the approximate version used below without 7.^^'' and self-energy corrections, the right hand side of the flow 
equation for the interaction is just given by the T-derivatives of the one-loop diagrams. 
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FIG. 2: a) The different fRG approaches in the parameter space spanned by interaction strength U , temperature T and 
infrared cutoff A. The surface represents the criticai manifold below which perturbation theory diverges, b) Box discretization 
of the BZ and the Fermi surface, c) A^'-patch discretization. The patches extend from the origin across the FS to the corners 
of the BZ. 



c. Interaction flow: In the interaction flow schemeQ we first multiply Q with a scale factor 1/ g and split it in 
two, yielding 



/2 



-iujn + e{k) ipkg 



(14) 



g will be the flow parameter. We can absorb the factor 1/g in rescaled flelds -0, ip defined as t/j 
the interaction term gets an extra factor g^ when written in terms of the new fields: 

1 



g-i/^V). With this 



(15) 



k,k' ,k + q 



We observe that changing the scale factor 1/g in Qg corresponds to changing the strength of the bare interactions. 

The rescaled fermions tp, 4> describe a system with a bare interaction strength g^ V . Now we can start aX g — 0+, i.e. 
at infinitesimally small bare interaction, and use the flow equations to integrate up to the desired bare interaction, 
reached at 5 = 1. We can also stop the flow at any other value of g, with the functions 5E and g'^Vg{k, k' , k + q) being 
the self-energy and interacting vertex function for the bare interaction g'^Vjk, k' ,k + q). We call this the interaction 
flow (IF). Related schemes have been proposed by PolonyiT^ and Medenp^. 

In the IF, singularities on the right hand side of the flow equation are not regularized by the flow parameter. Thus, 
the IF scheme has to be performed at T > 0, when the individual one-loop diagrams are bounded. The strength of 
the bare interaction is increased continuously, and in the course of the interaction flow potential singularities are 
approached from below. 



Fig.|21a) illustrates how the various methods detect perturbative singularities from different directions in parameter 
space. We emphasize that the T-flow and the IF scheme do not correspond to viewing a system on different length 
scales. 



III. APPLICATION TO THE 2D HUBBARD MODEL 

Here we review results [13 m obtained for the 2D t-t' Hubbard model with next-nearest neighbor hopping i' 7^ 0. 
The case t' — was studied with related cutoff- RG schemes e.g. by Zanchi and Schulz 3] and Halboth and MetznerQ. 



A. Calculational issues and approximations 

All fRG approaches to the Hubbard model have so far neglected the frequency dependence of the interaction 
vertex by calculating the static part and approximating all other vertices with the static one. This can be argued 
to be good as the initial interactions are frequency independent, and the corrections to the leading divergences 
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FIG. 3: a) Critical temperatures for the flow to strong coupling and dominant correlations at the van Hove filling, obtained 
with T-flow RG for A'^ = 48 patches and U — 3t. b) Critical bare interaction strengths for the flow to strong coupling, obtained 
with IF scheme at T = O.OOlt and the van Hove filling and 12 x 12 box discretization. 

in the static components are small by power counting. The next approximation is the neglect of the self-energy 
corrections. The consistency of this can be checked by calculating the self-energy from the flow of the couplings. With 
these approximations and the initial truncation of the hierarchy one is left with the flow of the coupling function 
14(^1, ^2, fca). This function of three wavevectors can either be discretized directly on a grid in the Brillouin zone 
(see Fig. |21b)), or, more physically motivated, in the so-called Af-patch scheme, first used in this context by Zanchi 
and SchulzQ. This scheme takes advantage of the fact that for the standard Fermi liquid instabilities the wavevector 
dependence of the coupling function V{ki, k2, fcs) has a relevant part, which is the dependence of the processes between 
different points kf , fcl", k^ on the FS when the kf are varied around it. Furthermore there is an irrelevant part, which 
is the dependence in the directions orthogonal to the FS. Hence, one calculates the V{ki, fc2, ^a) for fci, ^2, ^3 on the 
FS and treats it as piecewise constant when ki,k2, and k^ move within elongated patches stretching from the origin 
of the BZ to the (±7r, ±7r)-points (see Fig. Elc)). 

Then, the RG flow is started at an initial scale Aq for the cutoff RG, at temperature Tq for the T-flow, or initial 
interaction strength g"^ — 0+ for the interaction flow, respectively. The initial interaction is chosen as Vg(fci , ^2) ^3) — U 
(= g^U for the IF). What is typically encountered at low T is a flow to strong coupling, where for a certain flow 
parameter £c one or several components of Ve{ki, ^2, ^3) become large. At that point the approximations break down, 
and the flow has to be stopped. Physical information is obtained from susceptibilities that can be calculated from 
the flow and directly from the fc-dcpendence of Vi{ki, k2, k^). For example, if only processes with ki + k2 = become 
large in magnitude (typically strongly attractive), the flow to strong coupling corresponds to a Cooper instability, 
indicating a dominant tendency to form a superconducting state in the parameter regime beyond the divergence. The 
symmetry of the superconducting pairing can be read off from the fci-dependence of the coupling function. 

B. Flow at van Hove filling 

First we consider the flow at the so-called van Hove (VH) flUing, where the FS is kept fixed at the VH singularities 
of the free dispersion at fc = (vr, 0) and (0, tt). We change the band filling by varying the next-nearest neighbor hopping 
t'. 

Using the T-flow scheme one finds three regimes ;'f^| . In the first regime for small \t'\ the FS is properly nested and 
there is a antiferromagnetic (AF) spin density wave regime. Then, for increasing \t'\, the d^2_j,2-wave superconducting 
tendencies get stronger and there is a regime where the corresponding susceptibility diverges most strongly. Around 
t' = — 0.33i, the flow to strong coupling is suppressed by orders of magnitude. On the other side of this potential 
quantum critical point the flow is dominated by ferromagnetic ordering tendencies. The same qualitative picture is 
found with the IF scheme0, although in this approach the d-wave tendencies seem to be weaker by a certain degree, 
most probably due to the unfavorable box-discretization (Fig. |21b)) used. Also, a parquet approach focusing on the 
VH regions gives the same order of instabilities T9J. On the contrary, the cutoff RG misses the ferromagnetic regime 
completely. The reason is the above-mentioned inability of the cutoff RG to treat small-(f particle-hole pairs on equal 
footing with other tendencies. Instead it produces a low- Tc Kohn-Luttinger instability in the singlet channel. 
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FIG. 4: Summary of the RG flow for fixed t' = — 0.3t obtained with A'^-patch cutoff RG. The left plot is for electron doping, 
i.e. band filling n > 1, while the right plot is for hole doping, n < 1. The band is half filled for /j, ~ —0.65t. The insets show 
the location of the Fermi surfaces for the different chemical potentials. Above the thick lines the flow to strong coupling is cut 
off by temperature. 

C. Flow with fixed t': saddle point regime 

Next we briefly sketch the tentative phase diagrams from the RG flow for fixed t' = — 0.3i as obtained with the 
cutoff RG[i1[3. These values of t' might be relevant for the high-Tc problem. However, the Mott insulating state of 
the high-Tc compounds can not be obtained with our weak coupling approach. Hence, our results cannot be directly 
compared to the cuprates. We can however read off some tendencies. For these parameter values there are still 
quantitative differences between the various RG schemes. These mainly concern the location of the boundaries and 
the numbers for critical scales. The order of the regimes however is the same for all approaches. 

First let us consider band fillings n > 1 which correspond to electron-doping in the high-Tc language. Here the 
RG finds just two phases when n is varied. One is an AF nesting regime with relatively large critical scales closer 
to half filling. Increasing the particle number reduces nesting, and finally (tt, tt) particle-hole become impossible at 
low T (see FS in the left plot of Fig. ^J. Then (or at somewhat lower doping, depending on the RG scheme) the AF 
tendencies become cut off. What remains is a low-scale d-wave superconducting instability. 

On the hole-doped side with n < I the FS is close to the VH points and the saddle point regions around (tt, 0) and 
(0, tt) play a dominant role. Between an AF nesting regime near half filling and a d-wave regime further away we 
find another regime, the so-called saddle point regime. Here several channels, most prominently the AF spin density 
wave, the d-wave Cooper, and also the d-density wave channel grow with comparable strength. All these channels 
reinforce each other mutually 16]. Hence, it is not obvious that the system goes into a long range ordered phase. A 
similar flow arises in Hubbard two-leg ladder systems. There the ground state is known to be a short-range correlated 
Mott insulator. Therefore it is conceivable that the saddle point regions of the 2D system will adopt a similar gapped 
strong coupling state with a truncated FS p^lT^ . We have elaborated on this idea recently by exact diagonalization 
of effective Hamiltonians suggested by the RG flow[l8j. 

IV. CONTINUATION OF THE FLOW THROUGH THE INSTABILITY 

So far we have only discussed RG flows where the effect of self-energy corrections has been neglected. This 
approximation will certainly break down when the couplings get large and it is desirable to include at least parts of 
the self-energy. 

There have been several fRG studies^ HI El HI El which gave some insights how the normal state self-energy 
is altered by the flow of the interactions. Yet in these works the RG flow still had to be stopped when the coupling 
constants became large. Physically one may expect that in many cases the opening of a gap in the fermionic spectrum 
should regularize this apparent divergence. Recently we have analyzed this process in the reduced BCS model using 
the IPl cutoff- RG 25]. In this model only particle pairs with zero total momentum interact with an attractive coupling 
constant. A theory with two- and four-point vertices only can be seen to be exact. The essential step is to include 
a small symmetry-breaking off-diagonal part Ao into the initial condition for the self-energy. Now also anomalous 
propagators appear in the flow when the cutoff is lowered. They generate a anomalous interaction vertex which has 
4 incoming or 4 outgoing lines. The flow equations for the self-energy and the interaction vertices as shown in Fig. ^ 
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can be integrated down to cutoff A = without a divergence of the interactions. A necessary point is a reorganization 
of the flow for the four-point vertex proposed recently by Kataninp^l- The A = 0-result for the anomalous self-energy 
approaches the BCS result for Aq 0. With initial Aq > the flow of the four-point vertex remains flnite. Only for 
Ao — > 0, a combination of normal and anomalous vertices diverges like I/Aq. This corresponds to the vanishing mass 
of the Goldstone boson when the U(l)-symmetry is restored. 

So far this scheme has only been applied to the reduced BCS model. It will however also work for more general 
cases where the interplay of particle-particle and particle-hole processes drives the system towards a superconducting 
instability. If a nonzero Aq is included in the initial conditions, the flow down to A = will converge to a state with 
a gapped spectrum and a controlled magnitude of the interactions. The strength of this scheme is that no mean-field 
decoupling of the interactions is required. In particular, the attractive interaction which goes into the equation for 
the anomalous self-energy is generated automatically during the flow to all orders in the bare interaction. This should 
allow for a much more precise investigation of the phase diagrams of weakly interacting fcrmion systems. 
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